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Abstract. We demonstrate that a recently published quasiclassical M0ller type 
approach [Geyer and Rost 2002, J. Phys. B 35 1479] can be used to overcome 
the problem of autoionization, which arises in classical trajectory calculations for 
many electron targets. In this method the target is stabilized dynamically by a 
backward-forward propagation scheme. We illustrate this refocusing and present 
total cross sections for single and double ionization of helium by electron impact. 
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Classical trajectory descriptions of atomic collisions and ionization processes have 
a long history, dating back to the pioneering work of Abrines and Percival The 
method has never become a mainstream tool, but has been used over the years for a 
variety of collision systems. This so called CTMC (Classical Trajectory Monte-Carlo) 
method was originally formulated with macroscopic point particles, scaled down to the 
dimensions of a real hydrogen atom, but it can be derived as a discretized treatment 
of the system's Liouville equation, too (see, e.g. (D). In the hydrogen target, for 
which CTMC was formulated initially, the single electron orbits around the nucleus 
on a Kepler ellipse. If this concept is extended, the resulting many electron atom is 
highly unstable: the electrons exchange energy and finally one of them ends up in a 
very tightly bound orbit and all the others are kicked out of the atom. This process 
is called autoionization, as it does not need any external perturbation. Therefore it 
is practically impossible to use such a classical many electron atom as a target in a 
CTMC calculation — it dissolves on its own before the approaching projectile has any 
chance to interact with it. 

Various attempts have been made to stabilize classical many electron atoms, so 
that they could be used in ionization and excitation calculations. These attempts 
range from neglecting the target electrons' interaction completely in the independent 
electron model, through highly symmetric initial configurations, which autoionize 
slightly slower up to additional momentum dependent potentials to incorporate 
the uncertainty relation These ansatze are then usable in CTMC calculations, but 
they either describe a modified scattering system or lead to inconsistencies. 



Letter to the Editor 



2 



Motivated by these and other shortcomings of the standard CTMC method we 
recently proposed a quasiclassical description for particle impact ionization, which is 
derived as an approximation to the quantum mechanical description the M0ller 
formulation of the quantum scattering operator is translated into the Wigner phase 
space formalism and finally approximated by setting fi. — > 0. The approximation 
procedure itself is well established and the resulting method is technically very similar 
to CTMC, as in both treatments the cross sections are evaluated by propagating 
classical trajectories. But nevertheless there are two fundamental differences: (i) 
the phase space description of the target's initial state is derived free of ambiguity 
from the quantum mechanical wave function without the need to artificially introduce 
quantization recipes for many electron atoms and (ii) the M0ller form of the scattering 
operator translates into a classical backward-forward-backward propagation scheme, 
which stabilizes arbitrary initial distributions. With this ansatz we could calculate 
fully differential cross sections for electron impact ionization of hydrogen, which 
essentially reproduce the experimental results over a wide range of energies and 
geometries 1^. We will now demonstrate that this approximation can also deal with 
an autoionizing classical helium target. 

As mentioned above our quasiclassical approximation is derived as the h = limit 
of the Wigner formulation of the quantum scattering operator S in the M0ller form 
(for details, please see [|j and references therein): 

s = nln+ with n±= lim u\t)UQ{t) 

t — *=poo 

The propagators U{t) = cxp[— iiJt] and C/o(i) finally translate into solving Hamilton's 
equations of motion for each of the (multi dimensional) discretization points of the 
initial distribution 

p{t = 0) = p\x,p) = A/'^W„(5(x - Xn)5{p -Pn)- 

n 

The weights ?«„ are the values of the Ti = Q limit of the Wigner transform of 
the initial state wave function at the discretization point: Wn = w^{xn,Pn)- It can be 
shown that the Wigner transform is only one special case to select the initial conditions; 
by modifying the underlying correspondence rule nearly arbitrary translations between 
wave functions and phase space distributions can be constructed [Q, |^ . We will later 
use this freedom to calculate cross sections with a more simple initial state distribution. 

According to the M0ller scheme each trajectory first is propagated backward in 
time under the asymptotic initial Hamiltonian Hq, i.e., with the interaction between 
target and projectile switched off. If not denoted otherwise we will in the following 
use atomic units. Its form is then: 

The subscript p denotes the projectile whereas the target electrons are labeled with 1 
and 2. The nucleus is set to have an infinite mass. 

When projectile and target are separated far enough, denoted symbolically by 
t — —oo, the interaction is switched on and the trajectory evolves forward again under 
the full Hamiltonian H through the collision at t = and on, until the fragments are 
well separated again. Then the fragments are brought back from t = oo independently, 
i.e., with the asymptotic final Hq, to the initial time t = 0. If the initial state is 
unstable under the classical propagation, as is the case with a helium target, then it 
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Figure 1. One electron energies of the three electrons in the course of an example 
trajectory. The two target electrons (broken curves) each start at —1.94 a.u. The 
propagation starts backwards with (|l|) until t = —100 a.u. and is then reversed. 
The forward propagation is performed with the full Hamiltonian H. The energy 
of the projectile (solid line) is shifted from Ep = 2 keV to Ep = I Pi + /P2. For 
further explanations please see text. 

autoionizes already during the first backward propagation. When the projectile-target 
interaction is added at the turning point t = —00, it is negligible first; the forward 
propagation effectively undoes the autoionization and the projectile encounters the 
refocused target in nearly its initial state. The autoionization still takes place, but 
now it is shifted away to t < 0, where it has no influence on the actual collision 
dynamics, which takes place around t = 0. There is consequently no need any more to 
neglect some part of the interactions or to introduce additional stabilizing potentials. 

In a quantum treatment the first backward propagation only contributes a phase 
shift, as the initial state is an eigenstate of Hq. The cross sections remain unchanged, 
if it is neglected. In the classical approximation, though, the target is not stationary 
during the backward and the forward propagations; but if both are performed, most of 
the error due to the approximation cancels and the target is effectively stationary with 
respect to the central time t = 0: it is this point in time, where the initial conditions 
are set up, where the collision takes place and where finally the cross sections are 
extracted. 

The stabilizing effect of this M0ller type backward-forward-backward scheme is 

demonstrated in figure |l|: there the one-electron energies £"„ = ^ — n = p,l,2, 
of the projectile (solid line) and the two target electrons (broken lines) during one 
trajectory, i.e., the evolution of one single discretization point, are plotted against the 
propagation time. The interaction energies between the electrons are not included. 
The propagation starts at t = and first proceeds backwards, here up to i = —100 
a.u., under Hq (Q), i.e., with target and projectile independent of each other. Then the 
propagation reverses and the full H is used to propagate back to t = and through 
the scattering event. The plot ends at "-|-c»", which is t = 100 a.u. here. In the final 
backward propagation the electrons are independent of each other, i.e., their energies 
do not change any more. Hence, we do not show this part of the trajectory in figure 

The target electrons were started in a symmetric configuration for the trajectory 
of figure Within less than one period one of them is kicked into a large orbit. 
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from which it returns after about 50 a.u. and kicks out the other electron. After the 
propagation is reversed it is clearly seen that for most of the now following forward 
propagation the energies of the target electrons "rewind" the backward propagation, 
i.e., the broken lines are symmetric with respect to t = —100 a.u. In other cases (not 
shown here) one of the electrons is even kicked out of the atom to a positive energy and 
only comes back due to the reversed propagation. The projectile's energy oi Ep = 2 
keV = 73.5 a.u. is shifted in the plot by Ep — {I Pi +IP2), i.e., the projectile is plotted 
to start at the negative energy of the target. Ep is constant during the backward and 
most of the forward propagation. It only starts to change shortly before t = 0, when 
the interaction with the target electrons increases. Then the symmetry of the target 
electron energy trajectories against t = —100 a.u. is broken and all three electron- 
electron interactions together determine the dynamics of the actual ionization event. 
In this plot both electrons are lifted to positive energies, i.e., double ionization occurs. 

Another difficulty in the classical description arises from the fact that the final 
state is normally of well defined energy. Since the Hamiltonian is conservative only 
that part of the initial Wigner distribution contributes which lies on this energy 
shell. Consequently, the quantum distribution character of the initial state is lost. 
To overcome this problem we proposed in |^ to evaluate the cross sections in terms of 
the energy transfer, which, for hydrogen targets, is equivalent to looking at only the 
projectile's energy. 

With the helium target a trajectory contributes to double ionization, when the 
projectile's energy loss AEp is bigger than the sum of the ionization potentials IPi 
and IP2 of the target and when both the target electrons have gained at least half of 
the total binding energy: 

- AEp > I Pi + IP2 and AEi,AE2 > ^ (2) 
A contribution to single ionization is consequently defined by 

- AEp > IPi, AEi > — ^— and AE2 < — — -. (3) 

Of course, this test has to be performed with the target electrons' energy transfers 
AEi and AE2 swapped, too. 

To verify the conditions (||) and (^) for double and single ionization we need the 
initial and the corresponding final state wave functions. The initial state is the same 
in all cases. We are therefore, as in the experiment, able to extract all physically 
feasible cross sections from the same set of final values of the propagated trajectories. 

As a first test of the performance and consistency of our approach with the helium 
target we calculated the absolute total cross sections cr^^^ for single and cr^^-' for double 
ionization. They are compared to absolute measurements by Shah et al in figure 

i 

For this first calculation we have used a simple initial distribution obtained from 
a product wave function for the ground state of the helium target 

'>P{ri,r2)^ — cxp(-Zri)exp(-Zr2) (4) 

TT 

with effective nuclear charge Z = Each of the single electron wave functions 

is then translated into a phase space distribution by multiplying its densities in 
coordinate and momentum space Both the wave function and the resulting 
phase space distribution have a total energy oi E = —2.85 a.u., slightly less than the 
experimental value of —2.904 a.u. 
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Figure 2. Total cross sections cr'^' for single and cr'^' for double ionization: 
mparison of our results (solid lines) with the experimental data of Shah et al . 
(open and filled circles). 

The single ionization cross section a^^^ reproduces the measured data on the level 
of accuracy that is typical for a single electron CTMC calculation, see, e.g., Q: the 
maximum occurs at a lower energy and is slightly higher than the experiment, while 
the high energy behavior follows the classical 1/E decay . The explicit treatment 
of both target electrons and of all interactions can, of course, not reintroduce quantum 
effects like tunneling. On the other hand the accessible phase space volume is much 
bigger than with only one active electron and one might fear that the dynamics "strays 
away" from the reaction path of single ionization, completely distorting the cross 
section — which obviously does not happen. In fact our result is in good agreement 
with an nCTMC calculation by Schultz et al 

The more interesting result is, of course, the double ionization cross section 
a^^^ : figure || shows the first classical trajectory result ever, in which the dynamics 
according to the correct full Hamiltonian without any modifications or additions was 
solved - simultaneously both for cr^^^ and cr'^'. The double ionization cross section 
has two regions of different correspondence with the experiment: for high energies it 
decays as much faster than the experimental data. This suggests that in our 

calculation in the high energy regime both electrons are ionized independently, each 
contributing a classical factor of 1/E, and not in a sequential event, which should 
decay approximately as ct^^-' jl^. For impact energies below 250 eV, on the other 
hand, the experiment is reproduced remarkably well, both in shape and in magnitude. 
In that region no microscopic quantum mechanical explanation has been proposed yet. 
It is known, though, that right above the threshold the final state is defined only by 
the long range and long time dynamics of the outgoing electrons, which can be well 
described classically |4j. The good agreement between our classical result and 
the measured data shows that up to a total energy of about twice the total binding 
energy of the helium target the main reaction paths are the classical ones. 

In this letter we have demonstrated for the first time that electron impact 
ionization of a two electron atom, i.e., helium, can be calculated within a classical 
trajectory Monte-Carlo framework with the full, unmodified helium Hamiltonian. 
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This has been achieved with a quasiclassical M0ller formahsm: the propagation scheme 
refocuses and stabilizes the autoionizing target. The total cross sections, extracted 
from the energy transfer, compare well with the experiment to within the limitations 
of the classical approximation. 

The next; more demanding level of tests will be to compare the differential cross 
sections to experimental results and finally to understand the dynamics of double 
ionization in the low energy regime. 

This work was funded by the Israel Science Foundation. 
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